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Abstract 

A complete characterization of the structure of nuclei can be obtained by combining informa¬ 
tion arising from inelastic scattering, Coulomb excitation and 7 —decay, together with one- and 
two-particle transfer reactions. In this way it is possible to probe the single-particle and collective 
components of the nuclear many-body wavefunction resulting from their mutual coupling and diag- 
onalising the low-energy Hamiltonian. We address the question of how accurately such a description 
can account for experimental observations. It is concluded that renormalizing empirically and on 
equal footing bare single-particle and collective motion in terms of self-energy (mass) and vertex 
corrections (screening), as well as particle-hole and pairing interactions through particle-vibration 
coupling allows theory to provide an overall, quantitative account of the data. 
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Nuclear structure is both a mature M and a very active field of research mu , and 
time seems ripe to attempt a balance of our present, quantitative understanding of it. Here 
we take up an aspect of this challenge and try to answer to the question: how accurately 
can theory predict structure observables in terms of single-particle and collective degrees of 
freedom and of their couplings? 

In pursuing this quest one of two paths can be taken: 1) select one nuclear property, 
for example the single-particle spectrum, and study it throughout the mass table 2) 
select a target nucleus A which has been fully characterized through inelastic scattering and 
Coulomb excitation {A{a,a')A*), together with one- {A{d,p)A + l,A{p,d)A — 1) and two- 
((H -|- 2{p,t)A, A{p,t)A — 2) particle transfer processes and study the associated, complete 
nuclear structure information involving the island of nuclei A, A±1 and H±2 in terms of the 
corresponding absolute differential cross sections and decay transition probabilities. Here 
we have chosen the second way and selected the group of nuclei ns, 119 , 120 , 121 ,i 22 gj^ involved 
in the characterization of the spherical, superfiuid ^^°Sn target nucleus. 

Single-particle and collective vibrations constitute the basis states. The calculations 
are implemented in terms of a SLy4 effective interaction and a ni 4 (^S'o)(= Argonne 

pairing potential. HFB provides an embodiment of the quasiparticle spectrum while QRPA 
a realization of density = 2+, 3“, 4+, 5“) and spin (2=*=, 3^, 4=*", 5^) modes. Taking into 
account renormalisation processes (self-energy, vertex corrections, phonon renormalization 
and phonon exchange) in terms of the particle-vibration coupling (PVC) mechanism, the 
dressed particles as well as the induced pairing interaction were calculated (see [7]; 
see also [SHT^jh Adding to the bare interaction the total pairing interaction 

Up-f'f was determined. With these elements, the Nambu-Gor’kov (NG) equation was solved 
selfconsistently using Green’s function techniques [T6H20], and the parameters characterizing 
the renormalized quasiparticle states obtained. Within this framework, the quasiparticle 
energies are given by E,^ = (e^ — . The renormalised single-particle energy 

iy — ep = Zy[{ey — ep) + E™®"], is written in terms of the HP energy and of the even 
part of the normal self-energy, the quantity Zy providing a measure of the single-particle 
character of the orbital u. The state dependent pairing gap Ay = obeys the 

generalized gap equation 
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where = ul + and being the odd part of the 

normal self-energy, while is the cn—mass. The resnlting valnes of are shown in 
Fig.l. The contribntion of and to Ay are abont eqnal, density modes leading 
to attractive contribntions which are partially cancelled ont by spin modes, as expected 
from general transformation properties of the associated operators entering the particle- 
vibration conpling vertices [2TH23]. Theory (SLY4 -l-QRPA-l- (PVC) REN-I-NG) provides 
a qnantitative account of the experimental value 1.45 MeV). It is to be noted 

that in carrying out the above calculations use has been made of empirically renormalized 
collective modes. This is because SLy4 leads to little collective density vibrations (cf. Table 

l, where, for concreteness, the bare QRPA results characterizing the low-lying 2+ of ^^‘^Sn 
are collected [7], see also [21]), in keeping with the associated value of the effective mass 0.7 

m. In fact, collectivity is closely associated with a density of levels (~ m*) consistent with 
an effective mass m* = rrii^mk/m ~ m. This is achieved by coupling the two-quasiparticle 
QRPA SLy4 solutions to 4qp doorway states made out of a 2qp uncorrelated component and 
an empirically tuned QRPA collective mode [2S| (see Fig. 2 of [7], cf. also [2S]), an example 
of the fact that in a consistent PVC renormalised description of the nuclear structure, one 
has to treat (dress), on equal footing, all degrees of freedom (i.e. single-particle and collective 
modes). In this way, not only self-energy but also vertex corrections are consistently included 
(sum rules conserving processes), and thus the ’’bare” QRPA mode is properly clothed, 
bringing theory in overall agreement with experiment (see Table 1, second and third lines; 
for details cf. |27| as well as m)- 

To test how robust the results displayed in Fig. 1 are, we have recalculated Ay as a func¬ 
tion of the three parameters { rrik, /3j{J'^) } |28][29||25|. The results (|A/jii /2 — 
associated with the lowest quasiparticle state hii /2 are displayed in Figs. 2(a),(b) and (c). 
They provide evidence of the fact that a description based on the renormalisation of single¬ 
particle states and collective modes through PVC leads to a well funnelled nuclear structure 
landscape [30], displaying a global minimum for values of the set of parameters { } close 
to the empirical values: effective mass 0.7m, bare pairing interaction strength 

consistent with ni4(^5'o)(Go ~ 0.22 MeV, see [22]) and quadrupole deformation parameter 
(/^2)o ~ W2)exp ~ 0.13. 

Similar conclusions can be drawn from the study of the dependence on Up'’’”® and ^2 of the 
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quasiparticle spectrum associated with the valence orbitals (hii/ 2 , ^ 3 / 2 , Si/ 2 , ( 77/2 and ^ 5 / 2 ), 
of the splitting of the multiplet of states (hii /2 ® 2 +)i 5 / 2 -- 7 / 2 - and of the 7 — decay spectrum 
following Coulomb excitation, as can be seen from Figs. 2(d)-(f), Figs.2(g)-(h) and Fig.3, 
respectively. 

Because of the PVC mechanism, the different valence quasiparticle states undergo renor¬ 
malization and fragmentation, phenomena which can be specifically probed with one-particle 
transfer reactions. In Fig. 4(a) we display the absolute differential cross sections associated 
with the reaction calculated making use of the spectroscopic amplitudes 

associated with the strongest populated fragments of the valence orbitals hii/ 2 , ^ 3 / 2 , si /2 and 
^ 5/2 and of global optical parameters, in comparison with the experimental data [33]. The¬ 
ory provides an overall, quantitative , account of the experimental hndings. To be noted 
that the agreement found between the summed absolute differential cross sections associated 
with the almost degenerate state 3/2’'' and 11 / 2 “ (experimentally non resolvable [ 2 S|, while 
theoretically separated by 100 keV), results from a subtle incoherent combination of the 
I = 2, dain/dfl peak at 6cm ~ 20° and of that of the I = 4 one at 6cm = 47°. 

In discussing the ^^°Sn(p, reaction we concentrate on the ^ 5/2 orbital, the most 

theoretically challenging of all of the valence single-particle strength functions. This is 
because this state, being further away from the Fermi energy (e^g^^ = —11.3 MeV, ep ~ 
— 8 MeV) than the other four valence orbitals (see Table 2), is embedded in a denser set 
of doorway states (of type Si /2 ® 2+, ^ 3/2 ® 2+, gii 2 ® 2+, hii /2 <8 3“, etc.), as compared to 
the other ones. Consequently, it can undergo accidental degeneracy and thus conspicuous 
fragmentation. As seen from Table 3, although the calculated summed cross sections {a = 
6.15 mb) agree, within experimental errors, with observation (7.93 ± 2 mb), theory predicts 
an essentially uniform fragmentation of the strength over an energy interval of ~ 760 keV, 
while the data [32] is consistent with a concentration of the strength at an energy close to 
that of the lowest theoretical 5/2+ level (1090 keV). 

In keeping with the above scenario we have shifted the bare single-particle energy 
by 600 keV ((-11.3 -|- 0.6) MeV= - 10.7 MeV), amounting to a 6% change in the fc—mass 
(i.e. from 0.7m to 0.74m, Table 2 ), and recalculated all the quantities discussed above. 
Making use of the corresponding nuclear structure results and of global optical parameters, 
the absolute differential cross sections associated with the 5/2+ states populated in the 
reaction ^‘^^Sn{p,dY^^Sn and lying below 2 MeV have been calculated. They are displayed 
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in Fig. 4(b) in comparison with the experimental data. Theory provides now a quantitative 
account of the experimental hndings. In particular, of the fact that the strength function 
is dominated by a single peak. With the 600 keV shift, it is predicted at an energy of 1050 
keV carrying 4.4 mb and it is observed at 1090 keV with a cross section of 5.35 ±1.3 mb. 
The resulting overall agreement between theory and experiment is further conhrmed by Fig. 
4(c) where the absolute value of the one-particle transfer strength function associated with 
the population of 5/2’'' states predicted by the calculation is compared with experiment. As 
observed in Fig. 4(d), also this two-dimensional projection of the multidimensional nuclear 
structure landscape is funnelled, testifying to the physical robustness of the hndings. 

An alternative approach to the one discussed above which leads to almost identical hnd¬ 
ings regarding the ^ 5/2 fragmentation, can be obtained by treating the energy of the hve 
valence orbitals as parameters to be optimized selfconsistently within the framework of the 
full NG calculations, so as to best reproduce the quasiparticle spectrum (Opt. Table 2 ). 
The results can be expressed in terms of a state-dependent fc—mass {{mk)v/'rn) ~ 0.74 ± a, 
with a = 0.07. 

Making again use of the ehective occupation numbers resulting from the solution of the 
NG equation, the two-nucleon spectroscopic amplitudes of the reactions tY^^Sn{gs) 

and Sn{p,ty‘^^Sn{gs) have been calculated. With the help of these quantities and of 
global optical parameters, the absolute diherential cross sections have been calculated in 
second-order DWBA taking into account successive and simultaneous transfer, properly 
corrected from non-orthogonality contributions [3l]. They are displayed in Figs 1(b) and 
1(c) in comparison with the experimental hndings [H^llHU] . Theory reproduces the absolute 
diherential cross sections associated with the ground state transitions within experimental 
errors. The calculations have been repeated for diherent values of the strength of the PVG 
associated with the most important collective vibrational mode, namely the lowest 2~^ as well 
as for diherent strengths of the bare pairing interaction. While the dependence of (J 2 n{p, t) is 
very weak with f32 (not shown), it is conspicuous with An example of such dependence 
is displayed in the inset to Fig. 1(a). Again, this two-dimensional section of the nuclear 
structure landscape is of a well funnelled character. Within this context, it is noted that a 
measure of the reliability with which theory can describe the nuclear structure is provided 
by the relative dimensionless standard deviations arei (equal to e.g. a{Eqp)/(Eqp) in the case 
of the quasiparticle spectrum) associated with each of the different observables, and taken 
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at the minimum of the nuclear structure landscape, as shown in Table 4. 

We conclude that a theoretical description of nuclear structure based on single-particle 
(mean held with ruk ~ 0.7m) and collective motion (QRPA) and on their interweaving con¬ 
trolled by the particle-vibration coupling mechanism and leading to renormalization of both 
types of nuclear excitations through mass (self-energy) and screening (vertex) corrections 
and induced pairing, can provide an overall quantitative account of the nuclear structure 
representative of a mass zone (group of nuclei displaying homogeneous properties like e.g. 
sphericity and superhuidity, likely circumscribed by phase transition domains). Allowing for 
a weak state dependence of the /c-mass, determined by optimising the energy of the valence 
single-particle orbitals to reproduce the quasiparticle spectrum, the theoretical description 
of the nuclear structure probed in terms of direct reaction absolute differential cross sec¬ 
tions and based on renormalized single-particle and collective degrees of freedom, becomes 
accurate within a 10% error level. The PVC mechanism is found to play a central role in 
achieving this result. Within this context, we note that pairing in typical superfluid nuclei 
lying along the stability valley like Ii8.ii9.i20,i2i,i22gj^ ^ origin, in which and 

contribute essentially equally to the pairing gap. The above considerations and protocols 
are not only transferable to the remaining superfluid Sn-isotopes (not considered explicitely 
in the present case), but also applicable to the quantitative description of other spherical, 
superfluid nuclear mass zones. 
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nuj 2 + (MeV) 

B(E2 t) (e^ fm^) 

/?2 

QRPA (SLy4) 

1.5 

890 

0.06 

QRPA REN 

0.9 

2150 

0.14 

Exp. 

1.2 

2030 

0.13 


TABLE I. Energy, reduced E2 transition strength and corresponding deformation parameter (32 
associated with the low-lying 2'^ state of calculated according to QRPA and empirically 

renormalized QRPA as explained in the text, are compared to the experimental values m- 


Orbital 

(MeV) 

{mk)u/m 

SLy4 

Opt. 

<^ 5/2 

-11.3 

-10.7 

0.74 

97/2 

-10.1 

-10.5 

0.67 

Sl/2 

-9.0 

-7.9 

0.80 

d3/2 

-8.5 

-7.1 

0.83 

^11/2 

-7.1 

-7.45 

0.67 


0.74 (a = 0.07) 


TABLE 11. Energy of the valence orbitals associated with SLy4 (nik = 0.7m) and those obtained 
by optimizing (Opt.) the NG quasiparticle spectrum to the data. In the last column the results 
labeled Opt. are parametrized in terms of an effective, state dependent fc-mass. 


9 



























Exp. 



Th. 

e* (keV) 

a (mb) 

da/dO. (mb/sr) 


e, (keV) 

a (mb) 

da/dVL (mb/sr) 

921 

0.63 

0.75 


1150 

1.80 

2.3 

1090 

5.35 

7.0 


1290 

1.20 

1.7 

1354 

1.66 

2.3 


1710 

0.25 

0.32 

1562 

0.13 

0.16 


1910 

2.90 

4.0 

1730 

0.16 

0.18 






7.93 ± 2 

10.39 



6.15 

8.32 


TABLE III. The most prominent experimental (theoretical) fragments of the ^ 5/2 single-particle 
state populated in the ^^'^Sn(p,d) ^^®Sn (5/2+) reaction. The energies are listed in the first (fourth) 
column, while the absolute cross sections are given in the second (fifth), integrated within the range 
2^ < 9cm < 55°, and third (sixth), peak cross section, {9cM)max ~ 17°. The data are from [32] . 


Observables 

SLy4 

d ^/2 shifted 

Opt. levels 

A 

10 (0.7%) 

10 (0.7 %) 

50 (3.5 %) 

Eqp 

190 (19%) 

160 (16%) 

45 (4.5 %) 

Mult, splitt. 

50 (7%) 

70 (10%) 

59 (8.4 %) 

^ 5/2 strength (centr.) 

200 ( 20 %) 

40 (4%) 

40 (4%) 

d^l 2 strength (width) 

160 ( 20 %) 

75 (9.3%) 

8 ( 1 %) 

B{E2) 

1.4 (14%) 

1.34 (13%) 

1.43 (14%) 

cr2n{p,t) 

0.6 (3%) 

0.6 (3%) 

0.6 (3%) 


TABLE IV. Mean square deviation a between the experimental data and the theoretical values 
taken at the minimum of the corresponding functions displayed in Figs. 1(a) (inset),2(c,e,h),3(b) 
and 4(d) in keV for the pairing gap, quasiparticle energies, multiplet splitting, centroid and width 
of the 5/2+ low-lying single-particle strength distribution. In single-particle units Bgp for the 7 - 
decay (B(E2) transition probabilities) and in mb for cr 2 n{p,t). In brackets the ratio u/L, called 
(Trei in the text, between a and the experimental range L of the corresponding quantities: 1.4 MeV 
(A), 1 MeV (Eqp), 700 keV (mult, splitting), 1 MeV {d ^/2 centroid), 809 keV (=1730- 921) keV 
(^ 5/2 width), 10 Bsp (B(E2)), 2250 mb {a 2 n{p,t)), is given. 
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FIG. 1. (a) State-dependent pairing gaps for the five valence orbitals of ^^°Sn (see Table II, second 
column). The value of A associated with the HFB solution of ?;i 4 (^S'o)) is indicated by an arrow 
labeled The pairing gaps calculated making use of the empirically renormalised density 

modes are shown in terms of open circles joined by a dashed line, while the corresponding results 
obtained including also spin modes, and thus corresponding to Ajy are shown by the solid dots 
joined by a continuous curve. The contributions and A^”*^ are displayed in terms of solid 

triangles and solid squares joined by dotted and by dashed lines respectively. (b,c) Calculated two- 
particle transfer absolute differential cross sections associated with the reactions ^^*^Sn(p,t)^^^Sn 
(gs) and ^‘^‘^Sn{p,t)^‘^^Sn{gs) (continuous curves) in comparison with experimental data (solid 
dots) [35l|36]. In the inset of (a), the absolute value of the deviation of the integrated theoretical 
absolute cross section from the experimental value in the case of the second reaction is given as a 
function of the strength of the bare pairing interaction (cf. |29]i. 


11 





























FIG. 2. Absolute value of the difference between the experimental pairing gap and the theoretical 
value of Eq.(l)), calculated: (a) as a function of the effective mass associated with 

different Skyrme forces (solid curve) [2H], the dashed curve displaying the gap obtained using the 
fixed set of valence levels eu{Opt.) (Table II); (b) as a function of the ratio G/Gq [29]; (c) as a 
function of the ratio /32/(/?2)o associated with the lowest quadruple mode of ^^°Sn [25|. (d) (color 
online) The lowest quasiparticle energy values obtained from the full calculation as explained in the 
text referred to the energy of the 3/2~^ state, in comparison with the experimental data, (e) Mean 
square deviation between the experimental and theoretical levels shown in (d). (f) Mean square 
deviation between the experimental and theoretical energies of the five valence levels, as a function 
of the ratio G/Gq. (g) (color online) The experimental energies of the members of the /iii /2 ® 
multiplet are compared with the theoretical values, calculated as a function of the ratio /32/(/32)o- 
(h) Mean square deviation between the experimental and theoretical energies of the members of 
the /iii /2 <8) 2+ multiplet shown in (g). 
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FIG. 3. Experimental B(E2) values [M] in units of Bgp, of the quadrupole 7 — decay following 
Sn{a, Sn* Coulomb excitation, connecting the low-lying states of ^^®Sn (left). Also given 
are the theoretical values calculated making use of the results of the full renormalised calculation 
as explained in the text (right). The energies are in MeV. Mean square deviation between the 
experimental transition strengths associated with E2 decay from the 5/2^ levels, and the theoretical 
values calculated as a function of the /32 parameter is given in the inset (upper left). 
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FIG. 4. (a) (color online) Absolute finite range, full recoil DWBA theoretical differential cross 

sections associated with the low-lying fragments of the hii/ 2 ) c^ 3 / 2 ) ■Si /2 ^ 5/2 valence states 

most strongly populated in the reaction ^^°Sn(d,p)^^^Sn, calculated with the help of state of the 
art optical potentials and Vnp interaction (I.J. Thompson, private communication), making use 
of NG structure input, in comparison with the experimental data m- It is of notice that the 
single-particle orbit in the SLy4 mean field potential has been shifted towards ep by 0.6 
MeV (see text), (b) ^^°Sn(p,d)^^®Sn (5/2+) absolute experimental differential cross sections [32] . 
together with the DWBA fit used in the analysis of the data (right panel) in comparison with the 
DWBA calculations (left panel) carried out as mentioned in (a) . (c) Gomparison of the calculated 
strength function S^/ 2 {(T{^‘^^Sn{p, dy^^Sn{5/2~^)+a{^‘^^Sn{d,p)^‘^^Sn{5/2~^))/E with experimental 
data derived from one-neutron transfer reactions [321133] . The peaks have been folded together 
with a Gaussian function of variance 0.25 MeV. (d) The difference between the centroid (width) 
of the experimental and of the calculated d ^/2 strength S ^/2 is shown as a function of the ratio 
in terms of the solid (dashed) curve. 
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